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complex. For the larger multiple moment k, the magnitude of imaginary part increases with the 
increase of b, which means that presence of Lorentz violation makes Dirac field damps more rapidly. 
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to zero, which is quite a different from the asymptotic Dirac quasinormal modes without Lorentz 
violation. 
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I. INTRODUCTION 

Since Lorentz invariance was discovered, it has been great of importance in many fields of the fundamental 
physics, such as the Einstein's special relativity theory, particle physics and high energy physics. However, the 
recent development of unified gauge theories and the observation of high energy cosmic rays Q] P] Q imply that 
Lorentz symmetry is only an approximate symmetry of nature and may be spontaneously broken in the more 
fundamental physics defined in a higher scale of energy. Obviously, the existence of Lorentz violation will make 
a great influence on the fundamental physics and lead to many subjects need to be reconsidered. Therefore, a 



22 



9] 



lfj. It offers a 



great deal of effort has been attracted to study Lorentz violation in the different fields [4]- 

One of interesting theory models with Lorentz violation is the Standard Model Extension 8] 
consistent theoretical framework which includes the standard model and allows for small violations of Lorentz 
and CPT symmetry. This small spontaneous breaking of Lorentz symmetry maybe arise from the presence 
of nonzero vacuum expectation values for Lorentz tensors defined in an underlying theory. A straightforward 
method of implementing Lorentz violation in the curve spacetime is to imagine the existence of a tensor field 
with a non- vanishing expectation value and couple this tensor to gravity or matter fields. In gravitational 
theories, vierbein formalism is used widely because it can build a link between the covariant components T^ v ... 
of a tensor field in a coordinate basis and the corresponding covariant components T a b... of the tensor field in 
a local Lorentz frame. The link can be described by T^ u ... = e a „e „ • • -T a ;,..., where vierbein e a M is defined by 
9p,u = r labG a ^e a v . Moreover, it is also found that vierbein formalism can deal with local Lorentz transformation 
and diffcomorphisms, which are two basic types of spacetime transformations in gravitational theories. Bluhm 



and Kostelecky II] make use of this tool and find that any violation of diffeomorphism invariance via vacuum 
values of vectors or tensors breaks local Lorentz invariance. The converse is also true. For example, if a 
spacetime vector 6 M acquires a fixed vacuum expectation value (b^), which breaks diffeomorphisms, then the 
associated local vector b a as given by contraction with the inverse of vierbein also acquires a fixed vacuum 
expectation value (b a )- The presence of quantity (b a ) breaks local Lorentz symmetry. Therefore, it is natural 
for us to adopt vierbein formalism to study of Lorentz violation in the curve spacetime. 

Adopting vierbein formalism, the fermion partial action S^p in Standard Model Extension can be explicitly 
expressed as 

s^l^^r^-W), (i) 



where e M is the inverse of the vierbein e a . The symbols F a and M * are 



and 



M* = m + m 5 75 + a^ al a + b^ albl a + -H^ a e\a ab 



(2) 



(3) 



The first terms of Eqs.Q and ([3]) lead to the usual Lorentz invariant kinetic term and mass for the Dirac 
field. The parameters a^ L , b^, c liUl d^, e M , f lll g\^ v , H^ v are Lorentz violating coefficients which arise from 
nonzero vacuum expectation values of tensor quantities and comprehensive describe effects of Lorentz violation 
on the behavior of particles coupling to these tensor fields. All of coefficients can be constrained as the real 
numbers if the action ([1]) is hermitian. In generally, they are functions of position. The two terms involving 
the couplings a M and & M are CPT odd, which have been extensively studied in connection with Lorentz- and 
CPT-violating probing experiments including comparative studies of cyclotron frequencies of trapped-atoms 

n n 

17j, clock comparison tests [18| . spectroscopic comparison of hydrogen and antihydrogen 



anomalous magnetic moment 



13 ] , analysis of muon 
2l| . and so on. 



20l | , study of macroscopic samples of spin-polarized solids 
On the other hand, black hole is another interesting object in the modern fundamental physics, ft is believed 
widely that the study of black hole may lead to a deeper understanding of the relationship among the general 
relative theory, quantum mechanics, thermodynamics and statistics. This means that black hole physics play 
an important role in the fundamental physics. However, at present whether black holes exist in our universe 
or not is still unclear. A recent investigation shows that quasinormal modes can provide a direct way to 
identify black hole existence in our universe because that they carry the characteristic information of black 



holes 



241 ]. Moreover, it is also found that quasinormal modes have a close connection with the AdS/CFT 



correspondence 



25] 



2G] 
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29 L Thus, the study of quasinormal modes in 
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271 ] and the loop quantum gravity 
black hole spacetimes has become appealing in recent years 

Since both Lorentz violation and the quasinormal modes are hot topics in physics at present, it is natural to 
raise a question whether Lorentz violation affects the quasinormal modes of black holes. From the action ([1]), 
we can obtain that due to presence of Lorentz violating coefficients Dirac equation must be modified and then 
its quasinormal frequencies in the black hole spacetimes should be changed. However, to my best knowledge it 
is still an open question how Lorentz violation affects the properties of quasinormal modes in the background 
of black hole spacetimes. Obviously, different Lorentz violating coefficients (i.e. different types of breaking 
of Lorentz symmetries) have different effects on quasinormal modes. In this paper, we just consider Dirac 
equation with a modified term containing Lorentz violating coefficient b which arises from the nonzero vacuum 



expectation value of a special axial vector field 6 M [9([10( and calculate quasinormal modes of massless Dirac 
fields in the Schwarzschild black hole spacetime. Our result shows that in this case both fundamental and high 
overtones quasinormal frequencies depend on the Lorentz violating coefficient b. 

The organization of this paper is as follows. In Sec. II, we derive the equation of massless Dirac field coupled 
with the special axial vector field in the Schwarzschild black hole spacetime. In Sec. Ill, we evaluate the 
fundamental overtones quasinormal frequencies of the Dirac perturbational field by using the third-order WKB 



approximation 
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46] 



471 ] . In Sec. IV, we adopt the monodromy technique 



491 ] and study the high-damped 



quasinormal frequencies. The last section is devoted to a summary. 

II. THE DIRAC EQUATION WITH LORENTZ VIOLATION ASSOCIATED WITH AN 

AXIAL- VECTOR b M FIELD 

According to the variation of ip in the action (|TJ) , we find that the massless Dirac equation only containing 
the CPT and Lorentz covariance breaking kinetic term associated with an axial-vector b^ field in the curve 
spacetime can be expressed as 



[H a e a AI (a M + r M )-v/ 757 a ]* = o, 



(4) 



where 
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r 



a 1 
-a 1 



75 = i7°7 1 7 2 7 3 



I 

1 



(5) 



Since the Lorentz violation is very small, it is reasonable for us to assume that the axial-vector 6 M field does not 
change the background metric. For convenience, we take 6 M as a non-zero timelike vector (^-,0,0,0), where b 
is a constant. In the Schwarzschild spacetime, the vierbein can be defined as 

1 . „s 



1 



2M 



= , r, rsmf 



(6) 



Setting \& = (1 — — ) * (sin 6) = $ and substituting Eq. ([6]) into Eq. |4]) , the Dirac equation (J4J) can be simplified 



as 



I _ 2M dt 



, 2M „fl 1 . 7 9 7 3 d 
r or r r ou r sm a oip 



If we define a tortoise coordinate 



=757 



$ = 0. 



(7) 



and the ansatz 



r* = r + 2Mln( f) 



(8) 



(9) 



with 



■jm 



( JW Y r 1/2 \ 



—m y-m+1/2 



/ V 2j+2 r ; > 



V 



j+m+lym+1/2 
2j+2 I 



1 , 1 



(10) 



we find that the case for (+) and (— ) in the functions and can be put together and Eq.Q can be 
rewritten as 



_d_ ( G 

dr x \ F 



r 



2M 
r 



I 
-I 



u 
-w 



(11) 



It is very easy to find the decoupled equations for variables F and G can be expressed as 

d 2 F 

^ + (c 2 - V X )F = 0, 



d 2 G 
drl 



+ {lo 2 - V 2 )G = 0, 



(12) 
(13) 



with 



Vl , 2 = ± d -^ + w 2 



\k\Ai 



|jfe|A3 =p (r - 3M) - 26] 



6 2 ± 26(r - 2M) 



(14) 



where W 



2M 



\ and A = r[r — 2M). It is obvious that the potentials V\ and V% are related 
to the coefficient b of Lorentz violations, which means that the Dirac quasinormal modes should depend on 
Lorentz violation. If we set 6 = 0, Eas. (fT2^1 - (fT^) give the results of general Dirac fields in the Schwarzschild 



black hole spacetime 37|. Moreover, it is well known that the potentials V\ and V 2 possess the same spectra of 
quasinormal frequencies because that they are supersymmetric partners derived from the same superpotential 
W. In the following, we therefore just make use of Eq.(fT2"|) to evaluate the Dirac quasinormal frequencies and 
write Vi asV. 



III. THE FUNDAMENTAL DIRAC QUASINORMAL MODES WITH LORENTZ VIOLATING 

COEFFICIENT b u 



In order to study the relationship between quasinormal frequencies and the coefficient 6, we can take M = 1 . 
Fig. 1 shows the variation of the effective potential with the coefficient b of Lorentz violations for fixed k = 3. 
From this figure we can find that as b increases, the peak value of the potential barrier gets lower and the 
location of the peak (r = r p ) moves along the right. 



FIG. 1: Variation of the effective potential for the massless Dirac field (k — 3) with the coefficient b of Lorentz violations 



Let us now evaluate the fundamental quasinormal frequencies for the massless Dirac field with Lorentz 
violations by using the third-order WKB potential approximation, a numerical method devised by Schutz, 



Will and Iyer 
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4G] 



471 ]. Due to its considerable accuracy for lower- lying modes, this method has been used 
extensively in evaluating quasinormal frequencies of various black holes. In this approximate method, the 
formula for the complex quasinormal frequencies lu is 



lu 2 = [V + (-21/ ") 1 /2 A ] - i(n + i)(-2<) 1 /2(i + n), 



(15) 



where 



A 



n 



1 Vr 



(4)' 



(-20/2 8 I 



v 



(-2V ") \6912 I Vq 



-V (77+ 188a 2 ) 



1 

384 

1 

288 



V v 



(4)' 



(51 + 100a 2 ) 



1 



(4)' 



2304 \ y n 



(67 + 68a 2 



(5)' 



(19 + 28a 2 ) 



Vr 



(6) ' 



288 I Vq 



(5 + 4a 2 



(16) 



and 



a = n + 



d s V 



n is overtone number. 

Substituting the effective potential V into the formula above, we can obtain the quasinormal frequencies for 
the Dirac fields with Lorentz violations. The fundament modes frequencies for k = 1 ~ 5 are list in the Table 1. 
From Fig. 2, we find that for fixed k the real part almost decrease linearly with the increase with b. Moreover, 
Fig. 3 tells us that the relationship between the magnitude of the imaginary parts and b is more complex. For 
larger k, we find that they increase as b increases, which means that presence of Lorentz violation makes Dirac 
oscillation damps more rapidly. 
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b 


u (fc = 1) 


LO (fc = 2) 


LO (fc = 3) 


u (fc = 4) 


lo (k = 5) 




0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 
1.0 


0.176452-0.100109i 
0.166462-0.100956i 
0.157373-0.101275i 
0.149036-0.101037i 
0.141230-0.100325i 
0.133666-0.099621i 
0.126319-0.101000i 
0.122890-0.114215i 


0.378627-0.096542i 
0.367984-0.096909i 
0.357770-0.097175i 
0.347972-0.097342i 
0.338572-0.097408i 
0.329548-0.097380i 
0.320878-0.097266i 
0.312532-0.097083i 
0.304478-0.096864i 
0.296676-0.096668i 
0.289080-0.096612i 


0.573685-0.096320i 
0.562847-0.096596i 
0.552294-0.096829i 
0.542020-0.097021i 
0.532020-0.097174i 
0.522288-0.097288i 
0.512819-0.097363i 
0.503605-0.097402i 
0.494639-0.097404i 
0.485914-0.097372i 
0.477423-0.097310i 


0.767194-0.096276i 
0.756264-0.096489i 
0.745548-0.096681i 
0.735043-0.096851i 
0.724745-0.097000i 
0.714652-0.097128i 
0.704761-0.097234i 
0.695069-0.097319i 
0.685573-0.097384i 
0.676269-0.097428i 
0.667154-0.097453i 


0.960215-0.096256i 
0.949237-0.09643U 
0.938430-0.096592i 
0.927792-0.096739i 
0.917321-0.096873i 
0.907015-0.096994i 
0.896874-0.097101i 
0.886894-0.097194i 
0.877074-0.097275i 
0.867413-0.097342i 
0.857907-0.097397i 



TABLE I: The fundamental overtones (n — 0) quasinormal frequencies of Dirac field with Lorentz violating kinetic term 
associated with an axial-vector fo M in the Schwarzschild black hole spacetime for fc = 1 ~ 5. 



0.6 . 



0.2 0. 



FIG. 2: Variation of the real parts of Dirac quasinormal frequencies with the coefficient 6 of Lorentz violations. 



IV. THE HIGH-DAMPED DIRAC QUASINORMAL MODES WITH LORENTZ VIOLATING 

COEFFICIENT b u 



Motivated by Hod conjecture[28(, a great deal of effort has been devoted to the study of the high-damped 
quasinormal modes in the different black hole spacetimes because that Hod's conjecture suggests that there 
maybe exist a connection between the high-damped quasinormal frequencies and quantum gravity. In this 



section, we adopt to the monodromy method 48] 49] and investigate the high-damped Dirac quasinormal modes 
with Lorentz violations in the Schwarzschild black hole spacetime. Our purpose is to probe whether Lorentz 
violations affects Hod conjecture. As in Ref. 



|, after selecting the contour L as shown in Fig. 4, we can 
calculate the global monodromy around the contour L. In the neighborhood of the event horizon r = 2M, the 
effective potential V and the solution of Eq. (fl"2f can be approximated as 



Since the only singularity of <j>(r) or e~ lwz inside the contour occurs at the point r = 2M. After a full clockwise 



V (2M)« . 

round trip, <f>(r) acquires a phase e s , while e~ tuJZ acquires a phase e~~ . So the coefficient of e" 
in the asymptotic of Mr) must be multiplied by e 



(2M") 4 
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FIG. 3: Variation of the imaginary parts of Dirac quasinormal frequencies with the coefficient b of Lorentz violations. 
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FIG. 4: The complex r-plane and the contour L. The regions with the hachures denote the area Re(x) < 0. 



Moreover, we find that the behaviors of the tortoise coordinate r* and of the effective potential V in Eq. (|12 
near the singular point r = are 



V 



r 

b 2 - 4Mb 



4r 2 



(18) 



where j = 1 — m ■ As in Ref . 



according to the local monodromy around the singular point r = 2M, we find 
that the monodromy around the contour L must multiply the coefficient of e~ luz by a factor — (1 + 2cos7rj). 

Comparing the local and global monodromy, we can obtain directly the high-damped quasinormal frequencies 
formula for Dirac fields with Lorentz violations in the Schwarzschild black hole spacetime 



; 2 _ b 2 
(2M) 4 



[1 + 2costt(1 



2M' 



n — > oo. 



(19) 



Where k is the surface gravity constant of the event horizon of the black hole. Comparing with mass M of 



black hole, parameter b is very small, thus the term ( 2 m) 4 m Eq.(|19|) can be neglected. The frequency formula 
for high-damped quasinormal modes can be further simplified as 

w = Tffln|[l+2coS7r(l- JL)]\ -i2mrT H , n^oo, (20) 

where Tjj is the Hawking temperature of Schwarzschild black hole. It is obvious that Lorentz violating coeffi- 
cient b affects the high-damped Dirac quasinormal frequencies. 

In the higher dimensional (D > 4, D is the dimensions of spacetime ) Schwarzschild black hole spacetime, 
the axial vector field can be taken as .D-component form ( , 0, ■ ■ • , ■ • • ,0). Similarly, we find that in this 

D-l 

case the behaviors of the tortoise coordinate r* and of the effective potential V near the singular point r = 
are 



2(D-2)M' 
b 2 -2(D-2)Mb _ I -j 2 

V ~ r 2(D-2) 4 r 2 ' ' 

where j = 1 — tjjizwm • R- e P ea ting above operations, we can obtain that the high-damped Dirac quasinormal 
frequencies in the D dimensional Schwarzschild black hole spacetime satisfy 

w = Tglnll + 2cos7r[I - — ^ ]| - i2nirT H , n -> oo. (22) 

It is shown that the asymptotic frequency formula of Dirac quasinormal modes can be also affected by the 
coefficient b of Lorentz violation in the D dimensional Schwarzschild black hole spacetime. When b — » 0, we 
find that the real parts of high-damped Dirac quasinormal frequencies in both cases become zero, which agrees 



43 1 



44j | . According to Hod's idea [50(, one can obtain 



with the result of Dirac field without Lorentz violations 
that classical ringing frequencies with an asymptotically vanishing real part correspond to virtual quanta and 
the corresponding Dirac transitions in Lorentz invariance are quantum mechanically forbidden. However, from 
the formulas (|20[) and (|22[) . we find that Dirac quantum transition is allowable in the Lorentz violation Frame. 
It implies that the emergence of Lorentz violation may change the quantum property of Dirac field. 

V. SUMMARY 

Adopting the third-order WKB approximation and monodromy methods, we investigated the quasinormal 
modes of Dirac fields with Lorentz violating term associated with a special axial- vector & M in the Schwarzschild 
black hole spacetime. We find that the coefficient b of Lorentz violation affects Dirac quasionrmal frequencies. 
At fundamental overtone, the real part decreases linearly as the parameter b increases. But the variation of the 



10 

imaginary part with b becomes more complex. For the larger multiple moment k, the magnitude of imaginary 
part increases with the increase of b, which means that presence of Lorentz violation makes Dirac field damps 
more rapidly. Since the imaginary part of quasinormal frequencies for large multiple number k can be well 
approximated as 36j uii ~ ^J^qm ^ n 5^' ^ e P oss ^ e reason for our result is that the presence of the axial 
timelike vector field b^ with nonzero vacuum expectation values may lead to the decrease of Newton's gravity 
constant G, which is possible in the Lorentz violating theories. Moreover, Dirac also supports that Newton's 
gravity constant G decrease as increase of the age of universe. At high overtones, the real parts depend on 
Lorentz violating coefficient b. As b — > 0, we find the real part of high-damped Dirac quasinormal frequencies 
becomes zero, which consists with the result of Dirac field without Lorentz violations. Moreover, our result 
also shows that the emergence of Lorentz violation may change the quantum property of Dirac field. The 
effects of other Lorentz violating coefficients in the action (fTJ) on quasinormal modes of black holes need to be 
investigated in the future. 

Acknowledgments 

We thank the referee for his/her quite useful and helpful comments and suggestions, which help deepen our 
understanding of Lorentz violation and quasinormal modes. This work was partially supported by NNSF of 
China, Ministry of Education of China and Shanghai Education Commission. R. K. Su's work was partially 
supported by the National Basic Research Project of China. S. B. Chen's work was partially supported by 
the Hunan Provincial Natural Science Foundation of China under Grant No.05JJ40012 and Scientific Research 
Fund of Hunan Normal University under Grant No. 22040639. 



11 



[1] Zatsepin G T and Kuzmin V A 1966 JETP Lett. 4 78 

[2] Takeda M et al 1998 Phys. Rev. Lett. 81 1163 

[3] Coleman S R and Glashow S L 1999 Phys. Rev. D 59 116008 

[4] Arkani-Hamed N, Cheng H C, Luty M A and Muko-hyama S 2004 J HEP 0405 074 

[5] Csaki C, Erlich J and Grojean C 2001 Nucl. Phys. B 604 312 

[6] Cline J M and L. Valcarcel, JHEP 0403 032 (2004). 

[7] Libanov M V and Rubakov V A 2005 Phys. Rev. D 72 123503 

[8] Kostelecky V A 2004 Phys. Rev. D 69 105009 

[9] Colladay D and Kostelecky V A 1997 Phys. Rev. D 55 6760 
[10] Colladay D and Kostelecky V A 1998 Phys. Rev. D 58 116002 
[11] Bluhm R and Kostelecky V A 2005 Phys. Rev. D 71 065008 
[12] Kanno S and Soda J 2006 Phys. Rev. D74 063505 

[13] Alfaro J, Andrianov A A, Cambiaso M, Giacconi P and Soldati R 2006 Phys. Lett. B 639 586-590 
[14] Belich H, Costa-Soares T, Ferreira Jr M M, Helayel-Neto J A and Mouchereck F M O 2006 Phys. Rev. D74 065009 
[15] Henson J 2006 Macroscopic observables and Lorentz violation in discrete quantum gravity Preprint gr-qc/0604040 
[16] Martinez M R and Piran T 2006 JCAP0604 006 

[17] Bluhm R, Kostelecky V A and Russell N 1997 Phys. Rev. Lett. 79 1432 

[18] Bluhm R, Kostelecky V A, Lane C D and Russell N 2002 Phys. Rev. Lett. 88 090801 

[19] Kostelecky V A and Russell N 1999 Phys. Rev. Lett. 82 2254 

[20] Bluhm R, Kostelecky V A and Lane C D 2000 Phys. Rev. Lett. 84 1098 

[21] Bluhm R and Kostelecky V A 2000 Phys. Rev. Lett. 84 1381 

[22] Ferreira Jr M M and Mouchereck F M O 2006 Influence of Lorentz- and CPT- violating terms on the Dirac equation 

Preprint |hep-ph/0601018| 
[23] Chandrasekhar S and Detweiler S L 1975 Proc. R. Soc. London A 344 441 
[24] Regge T and Wheeler J A 1957 Phys. Rev. 108 1063 
[25] Witten E 1998 Adv. Theor. Math. Phys.2 253 
[26] Maldacena J 1998 Adv. Theor. Math. Phys.2 231 
[27] Kalyana S R and Sathiapalan B 1999 Mod Phys. Lett. A 14 2635 
[28] Hod S 1998 Phy. Rev. Lett. 81 4293 
[29] Dreyer O 2003 Phy. Rev. Lett. 90 081301 
[30] Cardoso V and Lemos J P S 2001 Phys. Rev. D 63 124015 
[31] Konoplya R A 2002 Phys. Rev. D 66 084007 
[32] Starinets A O 2002 Phys. Rev. D 66 124013 

[33] Setare M R 2003 Class. Quant. Grav. 21 1453 ; 2004 Phys. Rev. D 69 044016 

[34] Natario J and Schiappa R 2004 Adv. Theor. Math. Phys. 8 1001-1131 

[35] Leaver E W Proc. R. Soc. Lond. A 402 285 (1985) ; 1986 Phys. Rev. D 34 384 

[36] Ferrari V and Mashhoon B 1984 Phys. Rev. D30 295-304 

[37] Cho H T 2003 Phys. Rev. D 68 024003 

[38] Shao C G, Wang B, Abdalla E and Su R K 2005 Phys. Rev. D 71 044003 

[39] Wang B, Lin C Y and Molina C 2004 Phys. Rev. D 70 064025 

[40] Du D P, Wang B and Su R K 2004 Phys. Rev. D 70 064024 

[41] Zhidenko A 2004 Class. Quant. Grav. 21 273 

[42] Jing J L 2004 Phys. Rev. D 69 084009 

[43] Cho H T 2006 Phys. Rev. D 73 024019 

[44] Jing J L 2004 Phys. Rev. D 70065004; 2005 Phys. Rev. D 71 124011; 2005 Phys. Rev. D 71 124006 
[45] Schutz B F and Will C M 1985 Astrophys. J. Lett. Ed. 291 L33 
[46] Iyer S and Will C M 1987 Phys. Rev.D 35 3621 
[47] Iyer S 1987 Phys. Rev. D 35 3632 

[48] Motl L and Neitzke A 2003 Adv. Theor. Math. Phys. 7 307 
[49] Andersson N 1991 Class. Quantum Grav. 10 L61 
[50] Hod S and Keshet U 2006 Phys. Rev. D 73 024003 



